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Abstract 

We generalize a previous construction of a fermiophobic model to the case of more than one extra 
W and Z gauge bosons. We focus in particular on the existence of screening configurations and their 
implication on the gauge boson mass spectrum. One of these configurations allows for the existence 
of a set of relatively light new gauge bosons, without violation of the quite restrictive bounds coming 
from the pnc parameter. The links with Bess and degenerate Bess models are also discussed. Also 
the signal given here by this more traditional gauge extension of the SM could help to disentangle 
it from the towers of Kaluza-Klein states over W and Z gauge bosons in extra dimensions. 
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Introduction 



It has been argued several times in the literature that the Standard Model (SM) should be considered 
as a low-energy realization of a more fundamental theory, but the lack of evidence for the existence of 
new particles at LEP1 energies, apart from those predicted by the SM (with the exception of the SM 
Higgs boson), imposes severe restrictions on its extensions. Moreover, the precision in the measurement 
of the ratio between the strength of the neutral and charged current couplings, the so-called pnc 
parameter (NC stands for neutral current) ||, rules out models with a too large contribution to this 
quantity. 

At one loop in perturbation theory the contributions to /?nc grow at most quadratically in the masses 
of the scalar particles running in the loop. The quadratic terms can rapidly lead to a disagreement 
of the prediction with the experimental data |||. Therefore the phenomenon of "screening", i.e. the 
exact cancellation at one loop of the contributions to /?NCi which are quadratic in the mass of the 
scalar particles, becomes an important ingredient to be fulfilled by any viable alternative to the SM. 
In general the suppression of the new contributions could also occur numerically, via a fine-tuning of 
the parameter of the model, but we will not consider this possibility. We want instead to investigate 
in detail the occurrence of a screening phenomenon, i.e. the conditions to obtain an exact analytical 
cancellation of the terms quadratic in the masses of the scalar particles. Then we will analyse the 
implications of this requirement in the gauge boson mass spectrum. 

An extension of the SM, based on an SU{2)i x SU{2)r x U(l) gauge group, has been presented 
in a previous paper |Q. It has been called fermiophobic (FP) because of the absence of coupling 
between the fermions and the gauge particles of the extra SU{2)r group. The scalar sector of this 
model is the minimal one sufficient to give mass to the vector bosons and to introduce a non-trivial 
interaction among them: two doublets and one bidoublet. No triplet is needed owing to the absence 
of right-handed neutrinos in the model. Although, in general, SU(2) gauge extensions of the SM are 
characterized by large deviations in the pnc parameter, two sets of configurations of the parameter space 
have been found explicitly [|J|, which yield a screening of scalar quadratic contributions. Moreover, the 
model in 0], apart from being anomaly- free and without FCNC at tree level, passes succesfully all 
possible phenomenological checks allowing for the existence of relatively light new gauge bosons. A 
clear signature of this model would be the discovery of a relatively light extra Z' mainly coupled to the 
hypercharge. 

Motivated by the good properties shown by the fermiophobic model proposed in [|j in passing all 
phenomenological tests, we have studied a more general extension, based on an SU (2) i x SU{2) n x U (l)y 
(with n > 1) gauge group, focusing in particular on the existence of screening configurations and their 
implications on the gauge boson mass spectrum. 

Although the relation between custodial symmetry and screening is known |l[ ||] , we have preferred 
to investigate the general expression of pnc directly. We have not invoked a priori any extra global 
symmetry, but we have studied the conditions under which the screening phenomenon occurs and found 
the corresponding configurations in the parameter space of the model. We will see, in detail, how the 
requirements imposed by these configurations are so strong that they allow us to get the exact mass 
spectrum of the new gauge particles. The most remarkable consequence of this model is the possibility 
of having sets of relatively light gauge particles without contradiction with the experimental data 
regarding p NC . 

In a series of papers, Casalbuoni et al. [|| have constructed a gauge extension of the SM that 
includes, besides the standard 7, Z and gauge vector bosons, also two new triplets of spin-1 
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particles. They have called the model degenerate Bess. The interest of that model is mainly in its 
decoupling properties @. The latter originate from an extra global symmetry that the model has 
when the gauge couplings are turned off and that is also responsible for the degeneracy of the vector 
bosons masses. The non-linear realization of this model is able to reproduce the Higgs-less SM in the 
limit of infinite mass of the new vector bosons. A recently proposed linear realization of this model 
[^| has two interesting limits depending on the choice of the parameters of the scalar potential: in 
one case it coincides with its non-linear version (degenerate Bess) and, for a different choice, with the 
SM with a light Higgs field. Obviously this implies that at least at low energies the model shares the 
phenomenological success of the SM. We have studied the choices in the parameter space necessary to 
find a contact between our model, Bess and degenerate Bess. 

Similar extensions of the SM have already been proposed in the literature ||. Some of these 
extensions (with n = 1) appeared in the framework of composite [1C] or supersymmetric models [11, 12]. 



More recently, in the context of extra dimensions |13|, the tower of Kaluza-Klein states over the W 
and Z gauge bosons would resemble an infinite tower of and gauge bosons. The different 

type of contributions of these states to Ap and the e parameters as compared to the ones given here 
for a more traditional gauge extension of the SM, could help to disentangle both cases. 

The paper is organized in the following way. In the first section we present the Lagrangian of the 
model. In section 2 we give the gauge boson mass matrices in the interacting basis for the charged and 
neutral sectors where the photon is identified. We then discuss in section 3 the spectrum of scalars: 
Goldstones and Higgses. Section 4 is devoted to the computation of /?no In section 5 we describe 
the screening configurations and the corresponding mass spectrum for the gauge particles and we also 
compute the whole set of e parameters at tree level for the second configuration. We discuss in section 
6 the link with degenerate Bess and finally we comment on other phenomenological constraints and 
give our conclusions. 

We will adhere in this paper to most of the conventions used in j|. 

1 The model 



We consider a Lagrangian based on a gauge group G = SU{2)l x SU{2) n x U(1)y- The kinetic and 
self-interaction terms of the vector bosons are given, for each sub-group of G, by the usual field-strength 
tensor. We introduce a complex scalar doublet &j for each SU(2) group, and a complex bidoublet 
for every combination of two distinct SU(2) groups present in G (where j,k = 0,1,..., n label the 
corresponding groups). We need at least one complex doublet for each SU(2) group, to give mass 
to the new gauge vector bosons, and to break completely G down to U(l) em . We introduce also the 
bidoublets because they induce a non-zero mixing among the different SU{2) groups of the theory, also 
when these groups are characterized by sensibly different energy scales. As in the case of n = 1, triplets 
are not needed. 

This choice seems to be a good compromise between the idea of studying the screening in a general 
kind of scalar sector and the need of introducing the minimal extension sufficient to get an acceptable 
mass spectrum for all gauge and fermionic particles, and also to generate an interaction between the 
different SU(2) groups. 

The Lagrangian of the Yang- Mills sector is: 

1 " 3 1 
£ym = - t 12 G ) vafir ~ J B ^ B>1V i (!) 

* j=0 i=l 
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where 

Gjw = d^W^-d.W^ + gje^Wj^W^ i,l,m = 1,2,3; j = 0, 1, 2, . . . , n 
Bnv = d^B v — d u B^ (2) 

are the usual field-strength tensors, Wj is the gauge field of the j th SU (2) group, while corresponds 
to the U(l)y group. 

The quantum numbers of the scalar multiplets are indicated by a list of n+2 elements corresponding 
to the different groups: 

(0,l,2,...,n,Y) (3) 



and are: 



J 



1, 

2' 

3 k 



(1, . . . , 2 ,...,1,—-) for each doublet $j 



(1,..., 2 2 , 0) for each bidoublet V jk (4) 

j<k and j, k = 0, 1,2, ...,n, 

where 1 and the ellipses stand for a singlet behaviour with respect to the corresponding SU (2) group of 
the list and the 2 for a doublet. The values — \ and are the hypercharge of doublets and bidoublets, 
respectively. 

Starting from these assignments, and therefore from the transformation rules of each field under a 
gauge transformation, we write the covariant derivatives of the scalar multiplets. 

D»$, s (5) 
w,th* 3 = (I) 

is the covariant derivative of the complex doublet where g is the hypercharge coupling constant, gj 
is the coupling constant of the j th SU (2) group, and r are the Pauli matrices. 

D»* jk = d^ ]k -i gj T --W^ jk + tg k ^ jk T --W^ (6) 

with = j 

is the covariant derivative of the complex bidoublet tyj k , gj and g k being respectively the coupling 
constants of the j th and of the k th SU{2) groups. The ratio between a gauge coupling gj and go will 
be denoted by 

x 3 = 9 > (7) 
5o 

where also xo = 1 and x n+ \ = y = g/go is introduced to shorten the notation in the following. 

Starting from the definition of these derivatives, we explicitly write the kinetic and the mass terms 
of the Lagrangian of the scalar sector, leaving the scalar potential generically indicated with V(4>). All 
the fields are given in the interaction basis: 

n n 

£sca = £ {D^tf (D^j) + J2 Tr{(D^ jk )^(D^ jk )} - V(cf>). (8) 

3=0 3,k=0 j<k 
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We must introduce a gauge-fixing term, and we choose it in the usual way, in order to remove the 
mixing terms between Goldstone and gauge vector bosons: 



2C 



n 3 

EE 

j=0 i=l 



^ 2 



2C 



E 

3=0 



( )/^-^($;.t($ i)o -($t )Q ^ 



k=0 
2 



2 



£ TV ( ^ <* ifc >o - (*] fe )o T ^- fc ) + E Tr ( *w-5T<*lj>o - ^ot^ 
fe=i+i 

l 



(9) 



where 



(*i>o 



V2\ 



and (*jfe)c 



o 





m 



are the vacuum expectation values (VEVs) of the various doublets and bidoublets and 



(10) 



k = %fc-(*ife)o- 



We choose to work in the Landau gauge (£ = 0). With respect to fermions, since we assume this model 
to be fermiophobic as in the unique gauge-invariant coupling that can be constructed between 
fermions and scalars involves only the $o doublet, so the couplings are the same as in the SM. A first 
positive consequence is that one automatically avoids having FCNC problems at tree level, in contrast 
to the usual left-right models. Also, no mass term is introduced for the right-handed neutrino. 
Finally, the Lagrangian of our model can be obtained by combining these various terms, i.e. 



£ = -CyM + £sca + £gf. 



(12) 



2 Gauge boson mass matrices 



The kinetic term of the scalar Lagrangian in eq.(|8|) generates, after the spontaneous symmetry breaking, 
a mass term for the charged (Wj with j = 0, 1, • • • .n) and neutral gauge bosons (Zj with j = 0, 1, • • • , n). 

2.1 Charged sector 

The charged Wp boson mass matrix is an (n + 1) x (n + 1) symmetric matrix. In a compact way it 
can be written as 



M, 



Cij 



^Xi-i^-i (<%(^-i + 8(0) " °U ~ i ~ l)2«? _li - 1) 4 i - li - 1) ) with i < j (13) 



where 9(i) stands for the Heaviside function (with 6(0) = 1), u?- = v( + v^^ 2 [] and 



s (k) = 53 

i=0i/fe-l 



"ife-1- 



(14) 



1 wy is obviously symmetric under the interchange of its indices. 
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The mass eigenvalues can be obtained from eq.(p^) by performing an orthogonal transformation 



2d 

C 



RcM 2 c R T c 



(15) 

where Mq is the charged diagonal mass matrix and Rc is the rotation matrix that transforms the 
interacting charged Wf~ fields into the physical mass eigenstates W, 



p± 



W[± 



Rc 



(16) 



V Wit J 



This rotation matrix Rc involves mixing angles. The specific form of the rotation matrix Rc 

and the mass eigenvalues will be given in section 5 for the screening configurations. 



2.2 Neutral sector 

The neutral mass matrix in the interacting basis (Wq, Wf, Wf , • • • , W%, B) is an (n + 2) x (n + 2) matrix 
given by 

r'2 



Mi 



gl 



Xi—lXj—i 



(<%K ? _l + s{i)) - 6(j - i - l)u-_ lj _ 1 J with i<j 



L Nij ! 

where, to keep the notation as compact as possible, we have defined 

and v n+i = 0- 



(17) 



(18) 



The diagonalization of this mass matrix is done in two steps. First, we perform a finite rotation with 
a unitary matrix U to identify the photon field. This transforms the matrix of eq. (|l7|) in a block form 
with the first row and column filled with zeros corresponding to the zero mass eigenvalue of the photon. 
Afterwards, a second rotation Rn is performed to obtain the mass eigenstates of the Zi fields. 

Accordingly, the transformations from the interacting (Wq, Wf, W% , • • • , W%, B) to the mass eigen- 
states basis (^4, Zq , Zf, Z.f, . . . , Z^) are 

0\ 





A 






( w$ \ 




/ 


A 


\ 




f 


A 






( 1 • 


..0 ... 




z 






Wf 














z 










Zi 






wi 






zf 








Zi 














= u 




and 








= Rn 






where Rn = 





Rn 




Z n ~\ 






wl 






yP 








Z n -i 






V o 




V 


Z n 


) 




\ B ) 






yP 


) 




{ 


Z n 


) 







J 



The matrix Rn has the same structure as the matrix Rc- 

The finite matrix U has been built in two steps. In the first step, the first row has been constructed 
with the constraint of being orthogonal to the mass matrix of eq.(|l7]). In a second step, all the other 
rows have been constructed to be orthogonal to the first one. A different choice of the orthogonal rows 
simply implies a redefinition of the mixing angles of the matrix R^. The resulting matrix U is 



P(x fc _x) , 
ik = il /(— 1) — ' 



: i— 1 



/(fc-1) 

f(k — 2) 



6(k-i)6(i-2) 



P(x fc _ 1 )P(x i _ 2 ) 
f(i - 2)/(i - 3) 



(19) 
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where the functions P{xi) and f(i) are defined as 



P{xi) = — — and f(i) 



\ 



s=n+l 

E p ^sY 

s=i+l 



(20) 



Notice that f(n + 1) = 0. 

This is a generalization for arbitrary n of the U matrix introduced in [14] for n = 1. Applying the 
transformation matrix U to Mjy 

= UMlU T (21) 

we obtain the following matrix elements: 



M 



Nij 



9l 



0(i - 2)9(j - 2)1 if 5 l3 (v]_ 2 + s(j - 1)) - 9(j - i - l)u]_ 2l _ 2 



n+2 

E 

r=j 



'i-2 



E U sr-1 



+ 



f(i ~ 2) 



+ 



s=0 

i-2 



P(x, 



+ 



f(j - 2)' 



i-2) 



P(x, 



3-V 



2 " r _lj_2 



(22) 



E U r-lj 



f-2 



0(i-(i + l))(^_ 2 + S (j-l)) 



with i < j 



with 



2*3 



Xi-2Xj-2 

P(xo) 4 



f(j - 2)f(i - 2) 



f(j - 3)/(i - 3) 



1 



^j-a f(i - 2)f(i - 3)f(j - 2)f(j - 3) 



(23) 



The matrix elements of Mn have been denoted by Mn to recall that they are not the final result, 
because it is still necessary to make the second rotation on the Z fields in order to get the mass 
eigenvalues: 

(24) 



M 2 N d = R N MfjR A N . 



The complete rotation matrix is denoted by 



K = R N xU. 



(25) 



The precise form of the mass eigenvalues and the explicit dependence on the VEVs and coupling 
constants of the rotation matrix Rn or 1Z will be given in section 5 for the screening configurations. 



3 Scalars 



In this section we analyse the spectrum of scalar particles of the theory, Goldstones and Higgses. In the 
interacting basis of scalar fields we have 4(n + 1) degrees of freedom (d.o.f.) coming from the doublets 
together with 4n(n + 1) d.o.f. from the bidoublets. Half of them are charged, while the others are 
neutral. 
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3.1 Charged Goldstone bosons and Higgses 

The charged d.o.f. combine in 2(n + 1) charged Goldstone bosons Gf (i = 0, . . . ,n) and 2n(n + 1) 
states corresponding to the physical charged Higsses Hf (i = n+l,...,n(n + 2)). 

We will start by identifying the Goldstone bosons that will be "eaten" up to give masses to their 
corresponding gauge bosons. 

From the Lagrangian in eq.(||), after performing the rotations of the gauge fields from the interacting 
to the mass eigenstates, eq.(|l6|), we identify each Goldstone boson as the coefficient of the term linear 
in the corresponding physical gauge boson. 

Then we introduce a unitarity matrix A that gives the projection of the interacting fields onto the 
mass eigenstates 



± 



±(01) 

<Pi 

±(01) 
</?2 



<fi ijl) 



± (n-ln) 
± (n-ln) 



A 



Gf 



Gt 

H rv+1 
H-n+2 



(n+2) / 



c\ 



go(Rc)i lvo 
gi{Rc)i2Vi 



gn(Rc)in + lVn 
Ni 

Ni 

-go(Rc)ii4° ) +9i(Rch2v < ° 1) 

N, 



gj(Rc)i 3 +iVi -9i(Rc)n+iV2 
Ni 

-9j(Rc)i l) +9i(Rc)n+iv{ 3 Z 
Ni 



t T> \ (n — ln) , * (n — ln) 

gn-ljRclinVj -gn(Rc)in + lV 2 

, N > 

{ r> \ (n — ln) . / \ (n — ln) 

-g n -l(Rc)inV 2 +gn(Rc)in + lV 1 

Ni 



(26) 



where the normalization factor iVj is given by 



Ni = g 



\ 



E 

s=0 



n-1 



E 

3=0 l=j+l 



Xj(R c )ij+iv{ j ^ - xi(Rc)n+iv^ l) \ +vi<^v 2 



(27) 



The rotation matrix A can be split into two sets of columns A = {C\,C^}. The first set C\ (given 
in eq.(|26|)), with i running from 1 to n + 1, corresponds precisely to the decomposition of the charged 
Goldstone boson Gf_ 1 with mass mfi l in terms of the interacting basis. The second set with 
a = n + 2, . . . , (n + l) 2 is the corresponding projection of the mass eigenstate H^_ 1 whose mass 
is denoted by m^-i- We do not give here the explicit form of since we do not need it to find 
the screening configurations. In order to construct it is just necessary to find a set of states 
orthonormal among themselves and to the Goldstone bosons. For an explicit example in the case 
n = 1, see [[|]. Notice that the choice of this basis of orthonormal states leaves extra freedom to 
introduce n(n + l)(n(n + 1) — l)/2 angles. 

2 Notice that we use a different index notation for the Hf and their mass than in Our first physical charged Higgs 
is named H^ +1 instead of of ||, to be able to get a compact expression for pnc in the next section. 
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3.2 Neutral Goldstone bosons and Higgses 



The neutral sector is organized in a set of (n + 1) neutral Goldstone bosons and (2n + l)(n + 1) neutral 
physical Higgses. We assume here, as in that there are no other sources of CP-violation apart from 
the phase of the Kobayashi-Maskawa matrix. This hypothesis allows us to split the neutral scalars in a 
CP-odd set (n(n + 1) neutral physical Higgses Hf i = n + 1, . . . , n{n + 2) and n + 1 Goldstone bosons 
G° i = 0, . . . , n) and a CP -even set ((n + 1) neutral Higsses) without mixing between the two sets. 

In a similar way to the charged case, the rotation of the neutral gauge bosons from the interacting 
to the mass eigenstate basis in the Lagrangian of eq.Q, using eq . (]25|) , provides us with the definition 
of the neutral Goldstone bosons in terms of the interacting fields. 

Also here we introduce a matrix for the physical CP-odd states and the neutral would-be Goldstone 
bosons (called C) that gives the projection of the interacting fields onto the mass eigenstates basis: 



Imc^Q 
Ini(/)5 



Im0° 
Im</?2 (01) 



x O(n-ln) 

im(p 1 

T O(n-ln) 

lm<^ 2 



C 



G° 



n n+2 



n(n+2) / 



Ni 
Ni 

{gnKj + 1 n + l—gTZj + 1 n + 2)""' 
Ni 

(goTli+ii-gilZi+i 2)fi 01) 
Ni 

-(ffo^i+i i-gi~R i+1 2)^' 
N 



\ 



(01) 



(g£gj±i i+i-gi'R-i+i i+i)v[ 3 1) 
N 

-(gjlli+i j+ i-giTZ i+1 i +1 )v { 2 j l ^ 
N 



(9n-l^i + ln~9ii^i+ln+l)"i" 1 
N 

/ T> -r? \ — 1 n) 
— \9n-l /<-i+l n—9nl<-i + 1 n + 1 ) v 2 

N 



(28) 



where 



Ki 



P(Xj-l 



/(-I) +(RN) ^f(j-2 J s=2 
with i,j = 1, + 2 and the normalization factor is 



fti ~ 1) 



n+2 



Y^e(j- S )(R N ) 



PjXj-^PjXs-j) 

7( fl -2)/(a-3) 



(29) 



N = 50. 



\ 



k=n 
k=0 



n— 1 n 



yKi +ln+2 ) 2 vl + Y^ Yl (xjTZ i+lj+1 - xtTli+xi+xfu 2 ^. (30) 



j=0 



We can split the matrix into two sets of columns C = {C^,C^}, with the same conventions for the 
column indices as in the charged case. The column C£ given in eq.(28) corresponds to the neutral 
Goldstone boson with mass Tnf_i (i = l,...,n + l). The column index a runs from n + 2 to 
(n + 1) 2 and labels the set of physical CP-odd neutral Higgses fl^-i with mass m^-i- The construction 
of the columns C§ follows the same rules as in the charged case. Also in this case one has the freedom 
to define n(n + l)(n(n + 1) — l)/2 angles. 
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There is also a third matrix, which we have called B, that relates the real part of the neutral fields 
with a set of CP-even (n + l) 2 physical Higgses with mass mf_i (i = 1, • • • , (n + l) 2 )- We do not give 
an explicit form for this matrix since in the screening configurations it is always taken to be equal to A 
and C, except for the rows concerning the fields ^2 whose sign in matrix B is reversed with respect 
to the other matrices due to the conventions used in the bidoublets sector. This requirement simply 
fixes the rotation angles in the matrix B, and we always have enough freedom to choose them. The 
first column Cg corresponds to the SM Higgs-like with a mass m®. 



4 Computation of p^c 



In this section we give the explicit expression of the leading contribution (quadratic in the scalar masses) 
to A/9nc = Pnc — 1 fo r arbitrary n. 

In a renormalizable theory, at one loop, the contribution to ApNC can grow, by power counting, 
at most quadratically. As we have said before, we are interested in models for which there are nat- 
ural configurations of the VEVs and of the masses of the scalars such that the potentially large m? 
contributions to ApNC cancel and we have a screening phenomenon. 

The pnc parameter is defined following Q as the ratio between the neutral and the charged current 
couplings, and therefore the expression of Ap^Co evaluated at q 2 = 0, is: 



ApNC 



$V(0) E z (0) 



M 2 W 



Ml 



(31) 



where £y(g 2 ) (with V = Z,W) stands for the transverse part of the self-energy corrections to the 
propagator of the vector boson V. The contribution coming from £ 7 ^(0) is absent since we have 
checked that it is exactly zero by the orthogonality properties of matrix A. 

We use this definition of pnCi at g 2 = 0, since we are interested in corrections that are leading, uni- 
versal and independent of the value of q 2 . The kind of diagrams that can give a quadratic contribution, 
in the masses of the scalar fields, to eq.(|3l|) are shown in Fig.l. Their explicit expression is: 



z m,- 



t[mi,mj] 



16vr 2 



16vr 2 



(pUV + 1 — lnm 2 
(C uv + 3/2) (ml + m 



m„ 



In m? 



rrij In m 2 ' 



TTVX 



where Cuv = V e — 7 + Log4-7r and any symmetry factor or coupling constant has been factorized out. 

If one expands the Lagrangian in eq.@, keeping the vertices involving two gauge bosons (W, Z 
and A) and two scalars and those involving one gauge and two scalars, the result is 



ApNC 




+ {A s 2l+lB s i m+ i — A s ll+lB s 2 m+1 )Pk 



n n—X n 

22Aj+U+lBj + i m +lPj + 
3=0 j=0 k=j+l 

2 

f[mf, m^} 



(Asii+iBgi 



m+l 



A K 1M B, 



s2l+l£>s2m+l) Pj 
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[a) 



(b) 



Figure 1: Feynman diagrams contributing to /?nCj a ) z[mi] and b) t[rrii,mj]. The tadpole diagram 
(not included) coming from the scalar potential only renormalizes the VEVs and is not relevant for our 
computation Jl], [f| . 



n(n+2) 

+ £ 

Lm=0 



i n— 1 n r 2 

£ A'+l m+l-Pj + £ £ £ "^ s / '+ l( " s / m+lPj 

\j=0 3=0 k=j+l L /=l 

1 \ 2 



(»4 s 2/+lC s l m +l +A s il + iC s 2 



f[mf,m c m } 



AM 



■ n(n+2) 
Lf2 ^ £ 



2 Z,m=0 



i n— 1 n 2 

Vj=o i=o fc=j+i /=i 



i n— 1 n 

+ ( y^^j+lW^j+lm+l^f,- + £ £ (^sli+l^slm+l - C s2 /+lSs2m+l) ^jfe J /[^m; m f] 



where we have used that 
with 



t[mi,mj] = -z[mi] - z[mj] + /[mj, m 3 -] 



f[mi,mj 



1, 2 2\ mfm] 



-(m +mj)- 2 l J 2 Log[-^] 



JJ 16vr 2 

We have introduced, for convenience, the indices Sf (/ = 1,2): 



— m,A rri: 
l j j _ 



si = n + 2 



n 



j + fc 



s 2 = si + 1; 



(32) 
(33) 
(34) 

(35) 



■si labels the fields (f^ Re^ and Im^ ^ , while S2 does the same for 0^ ^ , Re</>2 and 



Im^ ° fc) . The functions P, JVf and AT± are defined as 

-fx = 5x (-Rc)lz+1, -Nj^. = 9jT^2j+l ± gT^2n+2, ^2jk = 9j^2j+l ± Qk^k+l- 



(36) 
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All the masses of the Goldstone bosons (mf- and mf, i = 0,1, ... ,n) should be taken equal to zero, 
since we are working in the Landau gauge. 

The expression of ApNC remarkably simplifies after the use of eq.(p3|) and of the ortogonality 



properties of the matrices A, B and C, independently of their explicit expression. From eq.(|32[) we 
notice that all the divergences cancel, indeed they cancel separately inside each self-energy; furthermore 
also all the finite contributions from z[rrii\ cancel. The function f[rrii,mj] gives a measure of the 
isospin breaking and its properties serve as a guideline to find the conditions to obtain the screening. 
In fact f[rrii,mi] = helps to find the conditions concerning the degeneracy of the scalar masses, 
while f[mi,0] = mf/(32?r 2 ) gives us the constraints on the VEVs in order to compensate the scalar 
contributions between the two self-energies. The particular conditions will be discussed in the next 
section. 

Also, we have verified that the previous expression reduces for the case n = 1 to the one of 0], in 
all the cases that are explicitly given there. 

From a direct inspection of eq. (|32|) , one notices the absence of terms proportional to f[mf, mf] and 



to f[m^ ,rrij], i.e. the absence of self-energy diagrams with two neutral scalars or two neutral physical 
pseudo-scalars running in the loop. On the other side, in the unitary gauge, the self-energy diagrams 
involving one physical scalar and one Goldstone boson are obviously absent and their contribution is 
shifted to the diagrams with the Goldstone boson replaced by the gauge field. From the previous two 
observations we learn that the contribution from the diagrams in Fig.l vanishes in the unitary gauge, 
if the scalar sector contains either only physical scalars or only physical pseudoscalars (cf. e.g. ||). 



5 Screening Configurations 



It was observed in j| for the case n = 1 that under certain conditions the contributions to /3nc> which 
are quadratic in the masses of the scalar particles, vanish. 

We have found that also when the gauge group is enlarged to SU(2)l X SU(2) n x U(l)y, there are 
at least two configurations in which these contributions to pnc vanish. The requirements imposed on 
the VEVs and consequently the predicted mass spectrum of the gauge bosons are one of the differences 
between the configurations. In the first one, which we will refer to as case I, the masses of all the new 
gauge bosons are sent to infinity, while in the second one, referred to as case II, a finite-mass spectrum 
for the new gauge particles is allowed. 

A common condition to both configurations is the requirement of alignment of the scalar transfor- 
mation matrices A, B and C or, to be more precise: 

A a/ 3 = B a/3 = C af3 except for B s2 p = -A s2 p = -C s2 p- (37) 

(jk) 

The difference in sign is just due to the convention chosen to write the 02 field inside the bidoublet 

ipjk- 

In order to understand the implications of this first condition on the parameter space, we should 
count the number of free parameters of the gauge and scalar sector (without the potential). We have 
on one side (n + l) 2 VEVs and n + 2 coupling constants and, on the other, n{n + l)(n(n + 1) — l)/2 
mixing angles among the charged scalars. In the neutral sector, we have n{n + l)(n(n + 1) — l)/2 
and n(n + 2)(n + l) 2 /2 mixing angles among the neutral pseudoscalar and the neutral scalar bosons, 
respectively. Of course, we also have the set of masses of the physical Higgses, already discussed. 
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This condition eq.(|37|) imposes constraints on the VEVs and requires each mixing angle of matrix 
A to be equal to the corresponding mixing angle in matrix B and C. So, at the end, in total, only 
n(n + l)(n(n + 1) — l)/2 free angles are left. Notice that the remaining mixing angles between the 
neutral scalars are completely determined by eq.([37|). 

A second condition will be required on the masses of the physical scalars. 

The particular requirements on the VEVs and on the scalar masses depend on the case and will be 
discussed in the next subsections. 



5.1 Case I 

In this first scenario, the set of equations given by (^) are satisfied by taking the large-Uj limit (with 
i > 0) while keeping vq, vf lj ^ and lj ^ fixed and small with respect to the Vi's. When a 
hierarchy between the v^s is introduced 

v < v i < v 2 < ... < v n (38) 

we find the expected alignment of the matrices A, B and C. In that limit the leading-Uj dependence of 
the charged and neutral rotation matrices is, respectively, 

(^o)ii~-j- and (■Rjv) ii ~^(i- 2)0(j" -2)4^i i<j, (39) 

where the elements along the diagonal are equal to 1 and the off-diagonal ones change sign when i > j. 
The leading term of the normalization factors is 



Ni ~ goJSa H + s(l)] +9(i-2)x 




2 2 
i-l v i-X 



Ni - 9o\l6a^ 7 ^H + s(l)]+9(i-2)xl 1 vl 1 ^-^. (40) 



Now it is easy to see, by looking at the charged matrix A and at eq. ([40]) , that in the large limit 
{k = 1, ...) the only non- vanishing term of the column i with i > 1 is the element in row i, which is 
equal to 1. All the other terms are suppressed either by a mixing angle or by the normalization factor. 
The first column i = 1 is exceptional since its vq is kept finite and the normalization factor does not 
grow. In that case all the terms proportional to a mixing angle are still vanishing while the others 
remain. Then this column turns out to be 



U A - 



vZ + s(l) 



/ n+2 row \ 

, ,o,...,o, ^,-4 oi) ,.. ! .r ) ,-4 on) ,o,...,o 



(41) 



V / 

It is not difficult to check, using eqs.(|28|) and (f40|), that the same holds exactly for the neutral sector. 
Notice that it is precisely the hierarchy of eq.(|38[) that guarantees that all the terms with a mixing 
angle vanish, since the VEVs in the numerator cannot be larger than the ones in the denominator. The 
standard gauge boson masses are in this case 



< f = f ^ + s ^> M k = 9 ijw^ + (42) 
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where we have used that /(0) 2 = /(— l) 2 — P(xq) 2 . The new gauge boson masses are, at first order in 
the large Vi expansion: 

2 

M k p = E i x i( v i+ s ( i+1 )) ( 43 ) 

„2 9l 2 /(^-l) 2 [ 2 lN PiPif ( 2 *v=? 2 //(i) 2 + /(i-l) 2 \ 2 

and go to infinity when the large Vi, i > 1 limit is taken. From the request that the coefficient in front 
of the electromagnetic current be the electric charge, one finds, once the matrix U has been applied on 
the interacting fields, that 

e = g P(x )/f(-l). (44) 



Using eq. (|44|) and the relation e = gsw, the sine and cosine of the weak mixing angle are respectively 
given by 

s w = P(x )/f(-l) and c w = /(0)//(-l). (45) 



Therefore, in this first configuration, it follows from eqs.fl42|) and ([4^) that p defined as 

p = M^/c 2 w M 2 z (46) 

is exactly 1. 

This configuration gives us the possibility to interpret in a simple way the alignment condition of 
eq, (|37|) . If one writes the Lagrangian of the model in this particular configuration, by using matrices 
^4, B and C, which now turn out to be equal and very simple, one sees that the Goldstones and physical 
scalar fields are organized in a transparent way: 

i=n i=n(n+l) 

C = Y J {D^i) ] D^i+ E (AA) t£ >/A (47) 
i=0 i=l 



where 



Ti = ( G i+_ ih ? \ and = ( H * + ^ h ]+n+l \ . (48) 



The set of Tj are complex doublets made of Goldstone bosons and a singlet Higgs, and fij are complex 
doublets of matter. Thanks to the alignment, Apnc is nothing else than the sum of the contributions 
of a set of complex doublets flq| . Therefore it, automatically, tells us which are the conditions on the 
scalar masses to obtain the screening configuration: 

raf" = mf or = mfi (49) 

where 1,1' = n + 1, ...,n(n + 2). We can now resume the main points of case I. In order to get the 
alignment condition eq,(^7|), we perform the large Vi (i > 0) limit, with a hierarchy among the VEVs, 
and we set equal each angle of matrix A to the corresponding of matrix B and C (so at the end only 
n(n + l)(n(n + 1) — l)/2 common angles are left free). Then we see that the only remaining source 
of scalar quadratic terms is given by the isospin breaking, i.e. the mass splitting in the doublets of 
eq. (|49^ ) . Once eq. (|49f) is fulfilled we have screening. 
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5.2 Case II 



This is the most interesting screening configuration. In fact a finite light mass spectrum is allowed for 
the gauge bosons. The result is obtained exactly, without the need of any approximation or limiting 
procedure. In this second scenario all the VEVs are kept finite. Obviously this makes it much more 
involved to find the solution of eq. (^7|) , so we will give some details of the way to proceed. We split the 
set of equations A = C into two sets. A first set is 



x k -x(Rc)ik Ni 

with k = 1, .., n + 1, while the second is 

XjTtj+ij+i — x{R-i+n+i (ji) _ Ni 

Xj vf\R c )i 3+ i - x lV f\Rc)u + i Vl N * 

XjH-i+ij+i - x{Ri + ii + i (ji) Ni 



(50) 



x j v^ l) (R c )i j +i-x l v < f l \Rc)ii+i 2 N i 
From the second set of equations one immediately obtains the condition 



(51) 



¥ t >=v$ ti . (52) 



Notice that our VEVs are taken to be real. Moreover, once this condition is imposed the fulfilment of 
eq. (|5l]) is automatically ensured by the fulfilment of eq. ( |50| ) . This is easily seen by making the difference 
between eq.(|50|), evaluated at k = j + 1, and the same equation evaluated at k = I + 1. It means that 
we can concentrate just on eq. (^0|) . We write eq.(|50|) in the following way: 



Xfc " l(jRc)ife **-i(*cr)ii with k<l (53) 



Xk-lTli+lk — yTli+ln+2 Xl-lR>i+ll — yTti+ln+2 

From the solution of the previous equation, we obtain the constraints on the VEVs, 

vi = — and ufj = uf i+1 ) , j > i. (54) 

Xi \ Xj ) 

The exact form of the rotation matrices (charged and neutral) corresponding to this screening 
configuration is the cornerstone of the calculation. The charged matrix becomes 



(Rc) ik = N, 



cl e(k-i) 



Xk-l 



(i - 2)5 ki -iXi_2 X] ~2 I > 



j=k x i 



(55) 



where 



n: 



c 



(56) 



l + d{i-2)xU E 



1 x. 
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while the neutral one is 



Hi k — Sn- 



P(x k - 



/(-l) 



(57) 



where 



with 



0(n + 1- k) 



Xk-l 



1 r 



8 n NF + e(i-2)9(k-i)Nf S ySii5kn+2N? + Xi- 2 Ski-iO(i-2)N t 



j=i-i x j 



rN 



N- 



N 



\ 



N, 



N 



N. 



c 



i > 2. 



(58) 



k=0 X k 



i + ir 



^ x 2 

j=0 X 3 



If we introduce these results into the general expression of ApNC displayed in eq.(^), the latter 
greatly simplifies since now Pj, and N 2 j k are constants, independent of j, k. It is then immediate 
to find the extra conditions necessary to get a vanishing A/9nc- 

a) The projection of the bidoublet fields on the physical set of Goldstones and Higgses should satisfy, 
for fixed I and any value of s\: 

(59) 



A 2 , , - A 2 



l S2 J+l> 

where I = n + 1, ...,n(n + 2) — 1, with the constraint that the sign of ^4 Sl z + i/M. S2 z + i should be the 
same for the whole column I. Notice that the constraint of eq.(^) is also satisfied in this configuration 
when I = 0, ...,n. A final important remark concerning this new constraint is that it corresponds to 
n(n + l)(n(n + 1) — l)/2 independent equations, exactly the same number of free scalar mixing angles 
as we have. Again, as happens for n = 1, this condition fixes the value of these angles. 

b) A second condition concerns, as in case I, the required degeneracy of the masses of the scalar 
particles in order to cancel A/?nc : 



A 



and 



m, 



m l 



c 



if A 



s 2 l+l/Asl l+l 



if A 



S 2 l+l 



i i+i 



+i 



l 



(60) 



(61) 



where I = n + 1, ...,n(n + 2). Notice that no restriction is imposed on the masses m® with j = 0, n 
(ttiq is the SM-Higgs-like). 

The conditions of case II are more restrictive than those of case I. In order to obtain the screening, 
here, we need not only to put constraints on the VEVs eq.(|5j) but also to fix the common mixing 
angles in the matrices A, B and C eq.(|5S), without any arbitrariness left. Also, as in the previous case, 
a constraint on the masses eq.(|6C|) and ( 6TJ) is required. 

Finally the spectrum of gauge boson masses turns out to be very simple in this second configuration, 



Ml P 



Ml? 



9l 



r It -■ 



+ 1 



9l 



^o 2 + E 

1=0 



2 2 i 2 
x l+l u ll+l + X i-1 



2 2 



-US. 



i > 0. 



(62) 
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It implies that, at tree level, p is not 1 but 



n 1 

A e = i-4Ei ( 63 ) 

5 = 1 X S 

where tw stands for the tangent of the weak mixing angle. As a consequence of eq.(|63|) it is clear that 
the model requires the extra gauge couplings to be large (as in [||) so as to obtain an acceptable value 
of p at tree level. However, large couplings do not necessarily imply large effects in the observables. It 
was noted in |4j that all the potentially large couplings in the charged sector are always suppressed by 
a small mixing angle (see eq.([55|) in the limit of large coupling with Xi > xj when i > j). On the other 
hand, in the neutral sector the terms proportional to a large coupling constant cancel in the sum of 
diagrams with light gauge bosons and heavy scalars together with the diagrams with only heavy gauge 
fields, leaving only subleading corrections to A/9no Similar scenarios were described in || and [15]. 



5.2.1 The e parameters: tree level 



Given the interesting spectrum of masses and properties of this second case, it is worthwhile to analyze 
it in more detail and compute the full set of e parameters ]i~8|] . We will follow closely the discussion in 
^ and the conventions for the parameters Ap, Arw and Ak. Similarly to what happen in [|| also our 
model exhibits a decoupling property in all e parameters. 

At tree level the Fermi constant Gf is given in our model in case II by 



Gf 
V2 



2<? 2 



E- 

s=0 



( 



\ 



2 



+ 



s=l X s 



J 



2vl 



+ o 



(64) 



where the first term in the large x\ expansion (l/2vg) would give the exact result if M w p = M w p. 

Following H we can also relate the cosinus of the weak mixing angle of our model with the corre- 
sponding one in the Standard Model denoted by cq in the limit of large 



4 



i-t 



i i 

1 1 32 + E ~Z2 

1 X s s= i x s 



(65) 



where 



2 1 
C »=2 + 



\ 



7TO 



V2G F Ml P 



(66) 



Comparing the SM definition of Arw 



Kf 



-Ar w 



(67) 



-w 
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with eq.(62) and eqs.(64-66), we find 

2 



(r 2 « 2 \ ( C VK S V/) ST 1 I 1 

— \ c w ~ s w) 2^Z2 3 ^2 1 

'W s=l « S W s =l x s y 



where we have used that M z p of eq. (|62| ) can be rewritten in a more convenient way 



(68) 



(69) 



In order to obtain Ak and Ap (current) we will compare with the standard definition of the Z current 



z o s c \ 2 



(70) 



where we denote with /j = {e^, e^, z/^, ujt, d^, <1r\ the chiral projections of the fermion fields and 
we leave implicit the generation indices. 

From the neutral-current (NC) gauge interactions of fermions 



w 3 
V B » J 



J Z () 



An 



( A, \ 



— ( Jem J\ 

one finds that the Z current in our model is given by 



( Jem 



( A, \ 

7 P 



V Z n» J 



(71) 



J Z P ~ ( R n)22 J Z +^2( R n)2+s2Jz s 



s=l 



that can be further decomposed using eq.(19) in 



4 



s 2 1 
bur 1 



CwJ k 3 Zy J * 

P(x n+ i)P(x s ) , L 
" f(s)f(s + l) B 



where 



(72) 



(73) 



4 3 =!?E 'l-u./r "/; and 3% = g £ Ffrffi 



(74) 
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Comparing eq. (f72|- 74 ) and eq.fl70|) one finally obtains 



1 + 



Ap 
2 

Ak 



S = l ' 



« =1 X « 



2+s2 



J x s f(s)f(s + l) 



E 



P(x ) 2 (^) 2+s2 



^ s 2 w x s f(s)f( S + 1) / - c w P(xq) 2 (^) 2+s2 ' 

^)22 + L Xamfi8 + 1) 

where (i?^-)2+ S 2 can be derived from eq.fli~9|) and eq.( 



(75) 



(R 



N)2+s2 



N 



N 



f(s-l)f(s)x s 



(76) 



In the large Xi expansion one finally gets 



Ap 



Ak 



1 A 1 



^2 ^ ™2 ^ ^ ^ ™2 I ^ 



c vy s =l X s s= i x s 



E 



i 

1 -^s 



(77) 



Recalling the definition for the e parameters in terms of the oblique corrections Ap, Ak and Aryy [18] 

ei = Ap 



c w ~ s w 



e 3 = c^/Ap + (c^ - s^/)A/c 
and using eq. (H) and eq.© we obtain 



(78) 



1 A 1 



«2 ^ ^ ™2 ^ ^ ^ ^ ff.1 I ^ 



1 n 1 f,2 _ 2 
r J- I l w * i. 



C2 



C3 



2 E „2 

s=l X « 



,2 



4+ I 1 



"3 + {c 2 w - s 2 w ) + 2 1 



(79) 



We find that similarly to what happens in the model described in [|| (and, of course, also in [||) there 
is a decoupling property in the whole set of e parameters in the large x^ limit. 



6 Links with Bess and degenerate Bess models 

It has been shown in Q that the FP model (n = 1) in the configuration II (v^ 1 ^ = v^ 1 ^ = v, v i = vq/xi) 
is able to reproduce the masses and mixings of the gauge boson particles of the Bess model [17] in a 
subset of the parameter space of both theories. 
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The translation table in our present notation is 



a -> 1/(1 + 2x1) 

f - vy/l + 2xl 

g" -> 250x1 

</ -> <7o2/- (80) 



In the left column, we have the 5 parameters of the Bess model (we use the notation of p7| ), which 
reduce to 4 free parameters with the extra constraint a = 2g 2 /(2g 2 + g 2 ). On the right, we have 
the 4 free parameters (go,x\,y,v) of our model after imposing the extra constraint v\ = \[2v. In this 
particular configuration the model shares all the good properties of the Bess model plus a screening of 
quadratic scalar contributions to Pno 

The link with the degenerate Bess Q, on the contrary, turns out to be more complicated, and we 
need to impose several constraints on both models. Beside the presence of an extra family of gauge 
bosons, the degenerate Bess model has also two important differences with respect to the Bess model: 

i) while the Bess does not decouple in the large mass limit, the degenerate Bess does; 

ii) the p parameter at tree level (eq.(^)) is different from 1 in the Bess model but equal 1 in the 
degenerate Bess. 

The degenerate Bess model has, concerning the gauge sector, 5 parameters (r, x, g, 9 or g' , v) in 
the notation of ||. Our model has on the contrary 13 free parameters in case I (although some of them 
become infinitely large), and 7 free parameters in case II. 

It is not difficult to see that case I cannot overlap with the degenerate Bess model from the condition 
that M\y 2 = My/ 1 . This is equivalent to imposing that x\v\ = x 2 v 2 , which is inconsistent in the large Vi 
limit with the necessary hierarchy found in case I. Even if one imposes an unnatural inverse hierarchy 
for the couplings, the rotation matrices of case I are no longer valid, breaking the alignment condition 



of eq.(37). On the contrary, case II fulfills this condition automatically. 

The imposed degeneracy between the new charged gauge bosons is translated into strong constraints 
in both models. In the degenerate Bess we need to impose 

x -» 0. (81) 

This condition implies a complete degeneracy of all new Z and W gauge bosons. On the other side, 
apart from the constraints of case II, i.e. 

V M = with ij = 01, 02, 12 

vo v 

V 2 = V\ = 

X 2 Xi 



5 2 

one should require in addition, in our model: 



«02 = «oA (82) 



«ia = x 2 > xi > 1, (83) 
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that is the corresponding of condition (^|). The second condition of eq.(^) is necessary to ensure the 

7 t 



absence of mixing between W^ 1 and the others W gauge bosons and also to get a tree-level p parameter 



equal to 1. 

The translation table of the remaining parameters of the degenerate Bess in the small-r limit is 

gv -> g v 
1 

ce -> , = c w 

v 1 + v 2 



vi + ui^i + xiy 



(84) 



where all terms of order 1/xf (i > 1) have been thrown away, and cyy stands for the large x% (i = 1, 2) 
limit of eq.fl45|) according to eq. (|83|) , The gauge boson masses are, according to eqs.(|62]), ( |8^ ) and 



M 2 - ^Ov, 2 M 2 - 9 ° v 2 

M W ~ d v M Z ~ d ~ 2 V 

MWi = M 2 W2 =M 2 Y =M 2 2 = 9 ^{vl + u 2 x (\ + x 2 )). (85) 
The new gauge boson masses are all degenerate and larger than the Mw mass in the small r limit. 



7 Discussion and Conclusions 

In this paper we have presented an extension of the SM based on a gauge group SU (2) l x SU (2) n x 
U(1)y (n > 1), with a non-trivial scalar sector. We have investigated in detail the /3nc parameter, 
with particular care for the scalar sector. Imposing the exact cancellation of the quadratic terms in 
the masses of the scalar fields leads to what we call a "screening" configuration. We have found two 
different possibilities for the screening to occur and we have calculated the respective mass spectra of 
the new gauge bosons. While in one case the screening phenomenon is obtained at the price of sending 
all the masses of the new gauge bosons to infinity, the second one is much more phenomenologically 
interesting, because a finite, relatively light mass spectrum is allowed, pnc severely constrains the 
parameter space of the model and we consider as a very good indication the fact that the screening 
occurs for arbitrary values of n. Indeed by analyzing the set of e parameters we have found that the 
model in case II exhibits a decoupling property in the large Xi limit. 

Of course, one should take into account other constraints coming from the precision measuraments 
at LEP and Fermilab. Indeed in a previous paper j3j, the phenomenology of this type of fermiophobic 
models was analyzed in detail in the case n = 1. A tree-level fit was done to the electroweak data, also 
the one- loop contributions to Rb, Ap and to all the flavour-changing processes that could be affected 
significantly by this model were computed. In that respect this kind of fermiophobic models (for any 
n) has a clear advantatge with respect to other gauge extensions due to the absence at tree level of 
FCNC and to the suppressed one-loop contributions to the process b — > S7 and to the Bq — Bq and 
Kq — Kq mixing. All these analysis show that relatively light W' and Z' are not excluded provided 
that x is sufficiently large. The Tevatron data could have changed this conclusion. However, analyzing 
the production cross section of the W 1 and Z 1 ', it was shown that the previous statement remains valid. 
Thanks to the fermiophobic nature of the new gauge bosons, the W couples to the fermions only 
through the charged mixing angle while the Z' has a direct coupling that scales as 1/x. Therefore the 
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weakness of these couplings that induce small cross sections allow the existence of relatively light new 
gauge bosons with no contradiction with the Tevatron data. 

All these conclusions can be applied also to our model, that reduces to the one studied in M in the 
case n = 1. One can see looking to eq.(|7^) that the contributions to the electroweak observables enter 
in an additive way. These contributions are organized in a tower structure, i.e. the lowest order W and 
Z families fix the bounds on the next ones. The inclusion of extra Ws and Z's receives then just more 
stringent bounds on the extra coupling constants Xi (with i > 1), mixing angles and masses. The good 
behaviour of the case n = 1 is not affected by the introduction of new gauge bosons; in the allowed 
region of the parameter space we can reach two different limit configurations: we can choose to have 
a very light new extra boson, letting all the others to be heavy, or to have a more balanced situation, 
with all the bosons in an intermediate mass range. 

The keypoint is indeed, as it was observed in J3J, the fermiophobic nature of these models that 
automatically allows to introduce a milder modification of the whole SM phenomenology, in contrast to 
other kind of gauge extensions of the SM. More precise and quantitative predictions can be obtained, as 
in the case n = 1 0], only through a complete phenomenological analysis, which is nevertheless beyond 
the scope of this paper. 
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